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STATE AND DEVELOPMENT OF FLUTTER CALCULATION 
Ey A. Telchmanti 


I. FUNDAMENTALS OF APPROACH TO THE PROBLEM 
A. General Treatment 


Schvarzmann.'s escposltlons eboved the flutter problem to be of 
such outstanding Importance that the large amount of vork for years 
expended on It seems more than Justified^ Nevertheless, the question 
arises of -whether It Is not sufficient to eliminate the danger of 
flut-ter by simply obeying a number of simple and generally known 
rules for the design, for Instance avoiding rearward positions of 
the center of gravl-ty. Identifying the separate natural frequencies 
of aerodynamlcally essential ccmponeht systems, obtaining rigid con- 
structions free from play, etc. 

The strict observance of such rules for design ca\xses, however. 

In many ways, structural limitations -which are by no means necessaiy. 
Figure 1 shows, for Instance, the critical speed of a -wing as a 
function of the natural frequency of the control surface. According 
to the rule of design, equality of frequencies bet-ween control surface 
and -wing (oi^ => cdq) ought absolutely to be avoided. Hb-wever, this 
-would be Justified only for the model at the right of figure 1, "model 
In the sense of the two-dimensional problem," -whereas for the "actual 
wing" at the left there Is no objection -whatsoever to equality of 
frequencies, since the critical speed there shows only an Insignificant 
reduction compared to the case of vanishing con-trol stiffness which 
must be taken Into account, anyway. 

On the other hand, ho-wever, there also exists the possibility of 
flutter occurring In spite of optimum observance of the design rules. 
Figure 2 shows the critical speed of a power-controlled auxlllary- 
con-troL-surface arrangement as a function of the position of the 
control- surface center of gravity. Even for mass-balanced control 
surface, and all the more for forward positions of the center of 
gravity, this system flutters at a relatively low critical velocity, 

■W'Stand und Entwlcklung der Flatterberechnung." Llllenthal-QesellBchaft 
f\ar Luftfahrtforschung Berlcht 135, PP- 11-20. (This paper -was 
presented at the conference on wing and tail-surface oscillations, 
March 6-8, I9IH.) 
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although the fln-hendlng frequency is 32 percent ahdve the auxiliary- 
control- surface frequency; and Just for rearward positions of the 
center of gravity the system exhibits a - insignificant, but still 
existing - region free from flutter. 

Examples of the type mentioned demonstrate the necessity of 
undertaJEing an Individually prepared investigation of its flutter 
characteristics for every airplane design. 


B. Individual Treatment 

In evaluating such individual investigations, it must be noted 
that for most parameters on which the flutter characteristics of an 
airplane design depend, reliable numerical values ere not available. 
This applies especially to control stiffnesses and control masses, 
to torsional stiffnesses of wing and fins which vary according to 
preloading and corresponding distortion (or wrinkling); above all, 
however, it applies, to the peirameters in the air-force law and to 
the structural -member damping. 

Moreover, the separate construction data for an alrpleme design 
change repeatedly, particularly at the initial steige where the flutter 
, investigation noist start if it is to affect the construction at all. 
Later on, too, many construction parameters fluctuate considerably 
due to tolerances or to subsequent modifications as they are required 
according to flight tests, that is, changes in the tall surfaces, 
lengthening of the fuselage, displacement of the control-surface 
center of gravity, variations of aerodynamic control-surface balance, 
etc. 


On the other hand, one must consider how abruptly a construction 
may pass from a flutter-safe to a flutter- dangerous state when a 
parameter is changed. (See fig. 3*) Accordingly, it would be funda- 
mentally wrong to "tailor" a flutter investigation only to one partic- 
ular co^lnation of individual design parameters without considering 
whether perhaps a parameter combination with distinctly unfavorable 
flutter characteristics is approached. A meaningful flutter inve^ti-_ 
gatlon must therefore comprehend, on principle, the entire reuiges” 
within which the uncertain design parameters may possibly lie. 

Of course, such "variational considerations" are necessary also 
when it is a question of finding the optimum design arrangement, for 
Instance the optimum position of split tall surfaces, the optimum 
degree of control- surface mass balance, or the most suitable eocrange- 
ment of an aerodynamic control-surface balEuace. 
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In view of the work expenditure on a conscientious flutter Investi- 
gation, researchers have prohahly pondered before whether It would not 
he advisable to set up detailed curve tables from which to take quickly 
the flutter characteristics of any construction. 

Unfortunately, there Is little prospect of realizing an "atlas of 
graphs" of that kind - even for quite special systems with two to four 
degrees of freedom the flutter characteristics depend on 10 to 36 freely 
disposable parameters. Besides, even If such atlases did exist, ex- 
tensive individual Investigations would still be necessary - merely 
the classification of a certain single construction within such an 
atlas would require the determination of decisive characteristics of 
that design, which Is an essential part of an Individual single In- 
vestigation. 

The uncertainties In the foiraulatlon for the air force seem to 
suggest, for flutter Investigations, that all calculations be omitted 
and all problems be Instead clarified by use of wind-tunnel models 
slmllaf to airplanes. However, one may call such tests superior to 
calculations only If the model has to a sufficient degree the aero- 
dynamic properties as well as the spring-mass and structural -member 
damping characteristics of the actual system. This requires, on one 
hand, performance of the teste at the correct Mach numbers, on the 
other, an Internal model structure which Is essentially simulated to 
the actual system and thus renders the latter ' s complicated bending, 
torsional, and camber characteristics - a modem airplane construction 
Is usually not replaceable by a stick with straight elastic axis. 

Also, the model must nowhere show higher s true turel -member damping 
coefficients than those corresponding to the actual system. In view 
of the technical difficulties In even approximately satisfying such 
requirements and guaranteeing their fulfillment, the construction of 
such a model progresses only slowly, according to all experiences. 
Consequently, the wind-tunnel test with correctly simulated models 
cannot be used for constructive decisions concerning the flutter 
problem at the design stage of an airplane, and such decisions can 
therefore be made In practice only according to calculations. 

Tests with correctly simulated models are intrinsically signifi- 
cant In that they make a conclusive checking of the flutter charac- 
teristics possible, but here again frequently extensive calculations 
are required, either because, alter all. In the end seme model data 
deviate from acttial conditions and thus require an additional calcu- 
lation for evaluating these deviations, or because unexpected results 
perhaps make a too high damping of the model or other discrepancies 
seem probable - aside from the fact that the Mach numbers of s\ich a 
test most probably will not agree with those for actual conditions. 

Of course, these objections do not refer to model tests of the funda- 
mental type, such as are necessary for obtaining data for detailed 
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calcTilatlonsj sucli ‘basic tests do not depend on the model's actually 
reproducing one particular airplane. 

The considerations above were meant to show from what viewpoints 
the DVL at present treats the .problems of flutter calculations. In 
closest collaboration with the airplane Industry, the attitude of the 
DVL proved to be right. Only In exceptional cases, single problems 
encountered by the Industry could be clarified Immediately from 
empirical facte, diagrams, etc; usually Individual calculations were 
required for the ultimate decision. 


II. NATURE OF IHE FLUTTER CALCULATION 
A. Mode of Consideration 


Because of the addition of the air forces to the Inertia, spring, 
and damping forces of the airplane structure, a flutter calculation 
must necessarily be more complicated than the otherwise Imown oscilla- 
tion calculations. Furthermore, modern airplane structures represent 
complicated combinations of framework and monocoque structures with 
cut-outs, hinged control surfaces, and auxiliary control surfaces, 
locally concentrated single masses, etc. Such strucfures can no 
longer be represented so as to be surveyable as easily as, for Instance 
sticks with straight elastic axis and well-defined bending and tor- 
sional stiffness. Accordingly, setting up differential or Integral 
equations of the flutter process generally does not by any means Imply 
a satisfactory answer to the technical problem. 

Since It Is Imperative to Introduce as few "unknowns" as possible 
Into technical, calculations. It would normally be Inexpedient to In- • 
dlvldually regard as unknown, for Instance, all elements of motion 
(paths, torsions, variations In camber) required In a complete descrip- 
tion of the flutter condition. Rather It will be advisable to con- 
centrate only on a few particular combinations Of all motion elements 
(compare fig. 4), that Is, certain "elementary forms" Fjj(x,z)j more 

will be said below on their selection (x,z = space coordinates). The 
xmknowns then would be the scales aj^ at which these elementary 

forniB F]j must be superimposed on one another. In order to describe 

approximately the actual flutter condition 6(x,z;t). At first, these 
scales aj^ are, of course, unknown functions of time. However 

it should be noted that in a motion according to forms F]^ prescribed 

in this manner a state of dynamic equilibrium can be established only 
when one takes Into consideration ihe additional forces which would 
have to attack at the system If it actually should carry out this form 
of motion. Hence, It Is advisable to use formulations stemming from 
the principle of virtual displacements for in case of ' displacements In 
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the sense of siiy forms the restraint forces do not produce vorh and thus 
drop out of the mathematical formulation. - The equations of motion for 
calculation of the time-dependent superlmposltlon scales aj^(t) of the 

separate elementary forme then may he directly set up fomlally In 

the form of Lagrange's equations. 

The difficulty In this treatment Is the selection of appropriate 
elementEoy forms. The natural oscillation modes of the system and 
modes somewhat similar to them are given preference In this respect. 

As Borhmann has shown recently, a group of n such modes F can he 
examined as to Its suitability hy determination of the air-force dis- 
tributions which would pertain to oscillations according to every 
single one of these modes F} the question then Is whether It Is 
possible. In turn, to approximate these air-force distributions by 
superlmposltlon of the assumed n modes F. 

Unfortronately, quite a nuiriber of such groups will usually be 
found appropriate according to the above procedure so that there Is 
practically no way other than to perform several flutter calculations 
with different choice of elementary modes until the miniimim critical 
velocity Is reliably determined. The lower natural oscillation inodes 
with few node lines are, of course, the most Interesting ones, since 
for them the Internal damping Is small compared to the absorbable 
air-force energy. 

In starting the design, one Is Interested first In wings and 
fins without control surfaces in order to make sure that the "basic 
structure" Is "all right" with regard to flutter. Next, the effects 
of control surfaces and auxiliary control surfaces, of controls and 
various fine details, are of Interest. 

The manlpiilatlon of flutter caJLculatlons would be made unneces- 
sarily still more complicated If one shoxild, for each separate calcu- 
lation, rigorously take Into consideration that the scale factors a]^(t) 

axe, at first, open functions of time. (Compare fig. ^. ) Hence, It Is 
advisable also to Interpret the Individual scale factors a]^(t) In 

turn as superlmposltlons of a series of pres.crlbed functions $]^(t), 

with g = 1 , 2 , . . . . (Compare fig. ^. ) The unknowns In the 
mathematical formulation then are, tilt 1 mat ely, the individual super- 
imposition scales of these time functions required 

determining equations are suitably set up formally with the use of 
Gauss' principle of minimum constraint. 

Concerning the question of which time functions are to be con- 
sidered for a technical flutter calculation, a general Indication 
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results, first, from the physically trivial fact that, for sufficiently 
small relative airspeed, every oscillation originally is necessarily 
damped, that therefore the transition^ to an excited oscillation must 
take place over a purely periodic intermediate state. It can he demon- 
strated that usually several such Intermediate states between damping 
and excitation or vice versa exist for a system. For the time being, 
only the lowest of these Intermediate states Is of technical Interest, 
thus a state which leads from the original damping to the first excita- 
tion.^ In practice, this Intermediate state Is mentioned as a harmonic 
state; that Is, cos cut and sin cut Etre selected as prescribed time 
functions 4jjg(t) and are suitably interpreted as components of the 

complex expresslcm It is, however, absolutely to be expected 

that In the future higher harmonic or any diminishing or Increasing 
time functions, respectively, will also be Included in the formulation. 

The advantages of a purely harmonic oscillation formulation csn 
of course be utilized only under the presupposition that harmonic forces 
pertain to harmonic motions. How far this Is admissible with respect 
to the actual air forces and to the stmctural member damping forces 
Is, at present, for want of systematic tests, an open question. How- 
ever, It certainly Involves considerable arbltxarlness In the case of 
systems containing a spring with Initial tension, that Is, of power- 
controlled control -surface mechanisms and of systems with free play. 


B. Flutter Equations 

After all these restrictions In posing the problem, the flutter 
equations may be set up formally as linear equation systems for the 
unknown superliiq)Osltlon scales a]^ of the separate elementary 

modes Fj^ (congwre fig. 5); then the 1^ equation, for instance, 

expresses the following: The Internal and external forces which 

would appear In flutter according to a linear combination of all 
elementary modes Fjj all together must not produce work If the 

system Is assumed deformed In the sense of the elementary mode Fj^. 
Correspondingly the separate coefficients of the flutter equa- 

tions are Integral expressions over products of the forces of a 
state ”k" and the motions of a state "1"; they axe thus analogous 
to the coefficients of the so-called elasticity equations, which 
are familiar to anybody working in statics research; they can be divided 
Into Inertia, spring, air, and damping members. 

therein the simplifying assumptions on which the air-force formulation 
Is based may cause certain difficulties; more details on this are 
contained In a research report by Lelss soon to be published. 
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Borljmanii recently vrote an exhaustive treatise on the nature of 
these Individual com.ponent terms (compare fig. and their determi- 
nation from calculations and tests. He therewith provided the long- 
needed continuation of various trains of thouf^t. Indicated hy the 
author In the Luf t f ahr tf or s chung 1939* A. Borkmann calls the s^arate 
components "mass, spring,' air force, and damping characteristics." 

Normally the mass siirface distribution will be known. The determi- 
nation of the mass characteristics then offers no difficulties whatsoever. 

Determination of the spring characteristics, however. Is less slmplej 
the quantity pjjj^ entering there (compare fig. 5) Is the loading which 
would have to attack on the airplane In order to force upon It exactly 
the deformation Determination of pj^^ Is simplest when the 

elementary mode F]^ Is precisely a natural oscillation mode of the 
system, because In that case equals the product of mass surface 

distribution, amplitude distribution F]^, and the square of the natural 
circular frequency pertaining to Fjj. 

In practice, however - above all In variational considerations - 
the elementaxy modes selected often will deviate i^re or less from the 
natural oscillation modes. Then It Is expedient, first, to determine 
separately the deformations pertaining to a series of load conditions, 
and then to superimpose the load conditions In such a manner that their 
resultant deformation agrees as well as possible with the elementary 
mode Fjj. The former superimposed loading then Is practically the 

required load.lng pj^ . 

Of course, the spring characteristics thus foiind are reliable only 
If In their determination the preloading existing In flutter which Is 
prodTiced by the steady lift of the airplane also Is taken Into considera- 
tion, for the preloading determines, wmong other things, the wrinkling 
and may thus cause considerable variations of the elastic properties. 

With the mass characteristics known, the spring characteristics 
may be determined also from the natural-oscillation conditions of the 
airplane. In practice, for Instance, on the basis of ground oscilla- 
tion tests which have to be performed. If necessary, with a preloading 
with soft springs. For this purpose, the Individual natural-oscillation 
mode first Is split up approximately Into components equal to the elemen- 
tary modes used. The resulting component coefficients of the Individual 
constituents are Inserted In the "static oscillation equations" as known 
amplltTides a. In order to then determine from them the spring charac- 
teristics contained In them. Correspondingly, one may of course deter- 
mine the mass characteristics If the spring characteristics have been 
determined otherwise. 
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Borkmaim, In hie aforementioned treatise, calls the methQds in- 
dicated "superimposing and splitting -up methods." All details con- 
cerning their nature and the manifold possibilities for their use 
are systematically treated in that report; of special interest are 
the poBslbllltles discussed therein for usln^ the methods for deter- 
mination of the air-force characteristics from oscillation tests in 
the air stream. I n as m uch as the mass and elastic characteristics 
have already been found in another way. Above all, iuch tests are 
significant in that they may be used for general statements on the 
air-force laws of the oscillating wing. Natureilly, simple models 
may be used to which the initially mentioned restrictions regarding 
the applicability of the model test do not apply. 

Flutter calculations of today are characterized by the linearized 
air-force formulations for the harmonicaJ.ly oscillating wing strip in 
the - strlpwise two-dimensional - flow of an ideal fluid such as were 
first developed by Kussner and later on eactended according to various 
viewpoints^ and adapted to practical use by Theodor sen, Ellenberger, 
Schwarz, Sohngen, and Dietze. The formulas at present ready for use 
deal, above all, with motions in broken and aeparated straight lines 
such as those sliown in figure 6. The imaginary models Indicated 
there represent various InteaTpretatlons regarding the effect of a 
wing strip with control surface and apply likewise to additional 
auxiliary control surfaces. Also it is not difficult to Insert, if 
necessary, curvatures into the form of motion. Dietze *s expositions® 
contain some ideas on the selection of the suitable imaginary model in 
the individual case. 

In calculations with these air-force formulations executed at 
the DVL, it proved desirable to be able to obtain all entering 
functions directly from curve tables. Thus the DVL recently in- 
stigated, for ma x imum time saving, the preparation of suitable curve 
tables uiLlform for all functions required in 1 interpretations of 
the wing with control surface and auxiliary control surface shown 
in figure 6. 

The few tests which, so far, give information on the validity of 
the customary air-force formulations show that for systems with control 
surfaces certain modifications are necessary even for smallest Mach 
numbers. Thus Voigt recommends, according to his flutter tests, 
certain reductions of the control-surface chord ratio used. 

It is quite certain by now that the customary air-force formula- 
tions no longer apply in cGLses of hig^r Mach numbers; therefore, modem 

®Dletze, F; Vergleichsrechnungen zum aerodynamlschen Ruderlnnenauaglelch. 

Lillenthal-Gesellschaft Berlcht 135, PP« 70-74. 
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flutter calculations should take the canpresslbility of the air Into 
consideration, all the more so since the formulations for Mach 
numbers =0.7 set up hy Posslo can he directly applied to arbitrary 
imaginary models. It is true that then for any imaginary model that 
happens to be in use, the numerical evaluation of a complicated 
integral, equation becomes necessary which expresses the adaptation 
of the flow variation to the dimensions and peculiarities of that 
imaginary model} this numerical evaluation must be made in succession 

for all reduced frequencies ^ ^ Mach numbers in question 

(cD = circular frequency, v = relative airspeed, I = wing chord) . 

Posslo has already carried out the evaluation to a larger extent for 
the wing. strip without control surface. 

Of course it would be valuable if one could so modify the foimiula- 
tlons customary since that time from the assumption of an ideal fluid 
that they can serve as manageable approximate formulae for congpreesible 
fluids} such attenqpta have been made occasionally, for instance with the 
use of Prandtl'e rule for steady flow. As far as is known, however, 
these endeavors failed or could not be generalized. 

It is also sure at present that the assumption of strix>wise two- 
dimensional flow may fall to work} according to tests and calculations 
by Cicala, for Instance, in case of so-called reduced frequencies 

= — — <il.O. Here also it is, for more accurate calculations, a 

question of numerically eval\iatlng an integral equation for every 
individual caae and for every required reduced frequency. Here also 
it will be correct to use first the formulations neglecting as little 
as possible and to turn to simplifications only after their admissi- 
bility has been proved. Accordingly, Kussner's and Posslo 's formula- 
tions are the first to be conslde 2 ;ed. Here, too, it would, of course, 
be valuable if one cotild so modify the expressions fpr stripwlse two- 
dimensional flow now in use that th^ wo\ild represent manageable approx! 
mate formulas for consideration of the induction. 

Independent of these analytical poBslbllitles for perfecting the 
air-force formulations for flutter calculations, it is, of course, 
imperative to check and develop them by means of tests. This includes, 
among other requirements, a careful clarification of their linear 
additivity. 

Thus it would be ideal if the Interference method developed by 
Zobel would soon be made disposable for air-force measurements on the 
oscillating wing. Independent of that, the DVL Intends to expedite 
indirect alr-rforce determination from model oscillation tests. 
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chetracterlzed previously In the discuasion of the conceptions "super- 
imposing and splitting-up methods." Voigt reports J on tests of that type. 

Far more uncertain than even the basis for the air-force charac- 
teristics are the data so f£Lr at disposal for the formulation of the 
damping characteristics. At the time it is not even. known whether 
it is at all permissible to linearly superimpose on one another the 
damping con 5 )onents corresponding to different motion components. 

Doubts concerning this fact arise from the obvious comparison of 
structuTEil -member dampings with friction forces which surely play a 
role in Itj for friction forces, a linear superimposition certainly 
is not possible. 

Thus, when today it is still regarded frequently as a rele, on 
principle, to neglect the structiiral-member damping in flutter ceLLcu- 
latlons, this fact may be Justifiable because this neglect will 
usually amount to a measure on the safe side unless the damping 
causes, exceptionally, a very far-reaching change in the coupling 
relationships. However, with the neglect of the structural- member 
damping, the actual nixmericaJ. representation of liie critical velocity 
is basically ruled out, and the flutter calculation merely permits 
separating the available designs into flutter-safe structures awfl 
into structures with flutter risk. Especially, recently there arose 
repeatedly the necessity of tnaki ng use of structural-member damping, 
possibly even artificial damping in order to Justify sensitive con- 
structions. Taking such cases into consideration, Boelk w-nfl Schmidt 
performed e^^rlments regarding the problem of damping years ago. 
Systematic investigation of structural-member flp.Tnp i pg often has been 
contemplated at the DVL. The first task will be to determine certain 
minimum values for the damping coefficients which in case of need may 
unhesitatingly be used in calculation. 


III. EATIOHALIZATION OP THE CALCULATION 
A. Slnqjllfled Air-Force Characteristics 


In view of the uncertainty of any determination of characteristics, 
the expenditure of flutter calculations is in an unfavorable proportion 
to their reliability. In calcxilation with the air-force formulations 
customary at present, a considerable part of the work is spent on the 
numerical evaluation of Integral e^cpresslons of the form 



^olgt, H.; Messung Instatlonarer Luftkrafte. 
Gesellschaft Berlcht 135, PP- 90-93- 


Lllienthal- 


V 
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(z = coordinate In direction of vldth, oi = circular velocity, v ■ rel- 
ative airspeed, Z ( z) « vlng chord at point z ) ; for these' expressions 

must he confuted f or nionerous values of the ratio (Con^are fig. J.) 


As Borhmann has shovn In the afore-mentioned report, this evaluation 
may unhesitatingly he shortened. If the transcendent functions appearing 
In the Integrand 


T 



are approximately replaced hy an • expression of the form 


ao + 

since then only the space function D(z) 'which Is Independent of os/v 
remains -under the Integral, sign. 

Leiss went e-ven further -when he r^laced — originally only for 
fundamental considerations — -the functions mentioned 


hy constant amounts. It -was shown that In this manner, too, practical 
flutter calculations can still he performed; however. It Is then advls- 
ahle for safety reasons to Intersperse occasional spot checks according 
to more accurate cadculatlons in-to e-very series of variations of the 
parameters of interest — In these spot checks Borkmann's simplifications 
may of co-urse he used. J'lgure 8 shows dependence curves calculated 
aocordlng to Leiss' simplification and Interspersed -wl-th numerous spot 
checks; -the figure gives an Idea of -the effect of Leiss* approximation. 
Naturally, only -the agreement at the heavily drawn sections of the 
lo-wer limiting curves of the flutter region Is of practlcsJ. Interest. 

Vlth Leiss' s Inpllf Icatlon and a spot check according to Borknann* s, 
the expenditure of a variational, consideration -with three degrees of 
freedom, for Instance, may he reduced from kj to I6 days (under the 
presupposition that all spring and mass characteristics are prescribed) . 



B. Method of Calculation 

The o-ther reason for the large -work expenditure of flutter calcu- 
lations lies In the search fot the critical -velocity after -fche matrix 
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of the flutter equations has been determined, the elements of idilch are 
known to contain the critical velocity and the flutter frequency cucr* 

The DVL was repeatedly faced with the necessliy of deciding on the 
fastest method of dete2milnlng the crltlceLl values v^ cu without taking 
Into consideration closely related secondary problems^ as for Instance 
the Increase of the flutter amplitudes after the crltlceLL velocity has 
been exceeded. 

On one hand, the critical state may be found from the flutter 
determinant (it then practically amounts to experimentation to find 
at which values of v, to the determinant disappears); on the other 
hand, "however, the critical state may be determined directly from the 
flutter equations; then It Is practically a matter of obtaining by 
Iteration the flutter form and therewith the critical velocity. 

The treatment of the determinant may consist In first splitting It 
up Into a real and a purely Imaginary conponent A ' and A ' ' and then 
plotting the amounts of the conponent determinants for different pairs 
of values of v, oi. 

For another treatment of the detesrmlnant, one may use the real 
notation. . _ 

Besides these primitive methods, special ones may be used, but onJy 
insofar as natural oscillation forma are selected as degrees of freedom, 
or as suitable transformations are performed first (which Is probably 
Isposslble at the outset In variational considerations of the type 
Initially described) ; furthermore, these methods require that the damping 
be neglected. 

We will first name a promising method, suggested by Borkmann in 
his afore-mentioned- report, as one of the conceivable possibilities of 
arriving at the pair of critical values v, co directly from the flutter 
equations; It consists In step-by-step Improvement of a component system 
selected from the flutter matrix by 6d.equate rise of all flutter equations. 
Another possibility Is to loprove an assumed flutter form ^by means of 
repeatedly carrying through the energy balance,^ 

In a report, as yet unfinished, Mayer has made estimates regarding 
the work expenditure for such methods. The points In figure 9 show a 
preliminary result; one is at liberty to double or halve the Indicated 
work periods. The curves shown -there are valid for calculations with 
formulation of the customary air-force laws and neglect, of the s-tiuctural 
member damping. One notices first of all -the great ad-vantage -the Iteration 
me-thods promise for more -than five degrees of freedom. Unfortunately 
this ad-vantage is coun-terbalanced by -the fact -that (due to -the entering 


^onpare Telchmann, Luftfahrtforschung (Aviation Research) 1939* 
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of the two parazneters v and <u) It la quite uncertain In the Individual 
case ■vdiether the method will converge tovard the group v, m of minimum 
velocity, or perhaps toward another group which la of no technical 
Interest, hut would then possibly he regarded as the decisive pair of 
values, For this reason, an entire region for the Iteration methods has 
been plotted In figure 9« Its lower limiting curves are valid for a cal- 
culation with a certain Initial system by carrying throu^ two Iteration 
processes; Its upper limiting curve shows what the expenditure would be 
If, as a precaution, three different Initial systems were tried out. 

lEie question of whether the minimum critical velocity of the systan 
under consideration Is actually found In this manner remains open, of 
course; besides, so far there exists no evidence whether and under what 
circumstances the Iteration methods for the flutter problem converge at 
all.. In view of the present state of knowledge one should forego the 
use of Iteration methods. If there are no reliable comparisons available. 

The more Involved methods which start directly from the flutter 
determinant are free from these disadvantages. In case -of less than six 
degrees of freedom. It Is there obviously advisable (according to fig. 9 ) 
to split the determinant Into Its real and Imaginary part; for more 
degree^ of freedom, however, the use of the real notation seems to be 
more advantageous with respect to time. 

As figure 10 shows, splitting up of the determinant Into a real 
and an Imaginary part Is advisable also In case of more than six degrees 
of freedom when Lelss' slnpllflcatlon of the air -force formulation (for 
Instance T' = T* ' = O) Is used. 

Figure 11 shows the work expenditure If the calculator uses Instead 
of the customary air-force formulations air forces determined, for Instance, 
by means of ezq>eriment so that varlotis theoretical relations do not enter. 

It Is Interesting that. In this case, use of the real determinant Is more 
expedient than use of the determinant split up Into real and Imaginary 
parte. 

At present. It Is still necessary to represent the combined effect 
of a larger number of degrees of freedom by separately considering numer- 
ous partial systems with two or three degrees of freedom each with the 
expectation that one of them svifflce to essentially describe the flutter 
characteristics of the entire system. In some cases, this might be more 
or less Justifiable by evalviatlon of the coupling relationship existing 
Just then. In many cases, however, some points remain uncleeo:; therefore, 
such a procedure should be \ised only as a last resort. This unfavorable 
situation can be changed only by considerably restricting the conputatlon 
times Indicated In the figures. In practice, nothing at all Is gained by 
suggesting that these times be reduced by 20 or even 50 percent; rather. 

It will be necessary to reduce them to' l/lO or l/20 of the amounts given 
In the figures. 
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This is poBslble only vhen au'toina'tlc calcula’blon devices are used 
such as are already employed for o"Uier purposes In America and England. 

In flutter problems such devices must accomplish two tasks; On one 
hand, they must detennlne the characteristics which requires above all 
the fonolng of multlt^xm product sumsj on the other h6uid, they must solve 
the flutter determinant or, respectively, the flutter equations for the 
desired pEtrameters v and Zuse now Is developing a device which Is 

to be suitable for both these tasks and many more. His device differs 
from those known so far by the fact that It e3<presses all numerical values 
In yes-no combinations of telephone relays; the course of the desired 
operations la then controlled by a prepared perforated tape. The D7L 
taken over the development of this device. 


DISCUSSION TO THE LBCTDHE OP A. TEICHMANN 


Bock .- It seems to me that one of the most Important viewpoints Is 
that the lecturer performs his calcttlatlons with consideration of different 
parameter variations and thus finds out Tdilch combinations of design 
parameters are unfavorable. For in the construction of alrplEine models 
the I n d i vidual parameters are not absolutely fixed but vary some^at 
from specimen to specimen so that avoidance of. such unfavorable parameter 
combinations Is lnq)ortant. 

Pvirthermore, I deem Important the efforts toward irestrlctlng the 
calctilation eaqjendlture, by mechanization of the calculation, so greatly 
that one may In the pTOject stage obtain really \isable surveys of the 
critical velocities as speedily as possible. 

flues sel. - Telchmann mentioned in his lecture, among other things, 
the sensitivity of the critical velocity with respect to the design 
parameters; therefore, I should like to state our opinion on this problem, 
and Illustrate It with a few examples. In particular. Just on the example 
given by Telchmann. We are dealing here not with laws of nature but with 
conclusions from a calculation experience Tdilch, of course, admits 
exceptions . 

A sensitivity of the critical velocity appears ordinarily only In 
the following cases; 

1. When in addition to a stubborn flutter possibility, a hsumless 
flutter possibility exists. 

By a stubborn flutter possibility, we understand one that is hard 
to el imin ate, by a harmless flutter possibility^ In contrast, one ■vdilch 
can easily be eliminated by relatively simple structureil measures. The 
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conception "hEirnileBS flutter posslbilily" has no connection with the 
conception of "henlgn" flutter. The harmlesB flutter posBlhlllty muBt 
he eliminated, ae a rule, because it cauBea too low critlceO. velocitiee. 
Then the Btuhhorn flutter posslhillty with higher valueB of critical 
velocity remainB which ordinarily 1b not BenBltlve with reBpect to the 
conBtruction parameterB. ThuB, a Bensitivlty no longer exlBtB after 
the required changeB concerning the flutter poBBihlllty of the airplane 
have been carried oiit. 

I Bhould like to take, aa an example, the caee named by Telcbmann 
for Illustration of the sensitivity; we Interpret It soimewhat differently. 
Figure 1(a) shows the variation of the critical velocity, calculated 
with three degrees of freedom bending- torsion-control surface, as a 
function of the control- surfeice imbalance with a seneitlvity of the 
critical velocity with respect to the aileron imbaJLance for rearward 
position a, as shown In Teichmann's lecture. (I sketched the figure 
from memory.) Figure 2(a) shows the diagram as it would appear 
according to our opinion If the dashed part of the curve were left 
out, the limitation to wiTilTmim values of critical velocity dropped. 
One recognizes that one deals here with two different types of flutter: 
first, the harmless bending^ flutter which, due to the agreement between 
flutter bending frequency, shows low values of critical velocity 
and a limited range of excitation and can easily be el im i na ted by 
damping or mass balance of the control svirface; second, the stiibbom 
or Intractable torsional flutter which remains after elimination of 
the bending flutter and shows, due to the approximate agreement 
between flutter and torsional frequency, higher values of critical 
velocity ftnfl no limited range of excitation. This second type of 
flutter may hardly be wholly eliminated. Is Insensitive with respect 
to design parameters, and therefore ordinarily permits only a moderate 
increase In critical velocity by structural measTires. 

2. When the critical velocity Is very high and, consequently, its 
dependence on a change In construction parameters Is not of Interest. 

We take as an exar^le the dependence of the critical velocity on 
the mass coupling for a wing with fixed ailerons (fig. 3 (a)) which 
shows a sensitivity of the critical velocity with respect to the mass 
coupling In the region a - b. Since flutter of three-dimensional 
wings with fixed aileron ordlnetrlly sets In - even for the most 
vinfavorable value of mass coupling - only at sufficiently high 
velocities, this sensitivity is not of Interest. 

3. When one deals with a physically Impossible calculation result 
which, due to the Imperfection of the assumed deformation lines, occurs 
in a region sensitive for calcTilation. 
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We take as an example the dependence of the critical velocity, 
calculated according to the power criterion, on the tending frequency 
for a three-dimensional wing with the degrees of freedom hendlng- 
tors ion- control surface (fig. ^(a)). In the region of curve 2, one 
could speak of a sensitivity of the critical velocity with respect 
to the tending frequency. As I es^lalned in detail in my contrltutlon 
to the discussion of Dlmpker's lecture, the critical velocity zero is 
Impossltle for a three-dimensional wing if even one single degree of 
freedom is elastic. In order to eliminate curve 2, one would have to 
perform an iteration, which is vei*y troutlesome . We introduced material 
damping. It wae found that a very STnall value of material damping was 
sufficient to eliminate curve 2. 

Thus a sensitivity of the critical velocity with respect to the 
construction parameters - the thing that matters - will hardly occtir. 

Telchmann expressed the opinion that it is not of Importance to 
calculate a single critical velocity hut to evaluate the safety 
against flutter within a larger scope tdilch is made possible, among 
other expedients, by variation of different construction petrameters. 
Furthermore, Telchmann pointed out that the number of parameters is 
vezy high. We are absolutely of the same opinion; therefore, we 
always carry out variations of the construction parameters. In 
selecting the parameters to be varied, we favor those with the 
following characteristics: 

1. The construction parameter has a great Influence on the 
critical velocity. 

2. The determination, of the parameter is unreliable. 

3. Variation of the parameter is convenient with respect to 
calculation technique. 

f 

In this manner, we arrive at a restricted number of variations 
in construction parameters and thus at a tolerable work expenditure. 

The intention of the DVL to follow up the problem of material 
damping is very commendable. For the time being, one is dependent 
on the tests compiled by Kussner. These tests result in a minimum 
value of the material damping a = 0.03. The number of airplanes 
investigated is not sufficiently large, most of them are not all- 
metal airplanes, and all of them eu:e old models. Although the ' 
material damping usually does not amount to much, knowledge of a 
reliable minimum value of material damping is very desirable for the 
following reasons: 
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1. If large znasseB participate In the motion, the Influence of 
Internal damping may hecome decisive. 

2. In certain frequency ranges the material damping may Increase 
the critical velocity veiy conaiderahly. 

3. Samettmes the calculation Eurrlves at physically Impossible 
results which disappear with Introduction of material damping, so 
that Involved Investigations become unnecessary. 

4. In exceptional cases, sensitivities occur when none of the 
presuppositions treated above are met. These sensitivities produce 
difficulties. As a rule, such, sensitivities cease with Introduction 
of a moderate material doping. 

Boelk . - Telchmaon mentioned a method of Borkmann according to 
which an Iterative solution of the determinant Is su gg ested, In such 
a manner that a partial system of the determinant Is solved with, 
respect to v,(U and this solution Is then Improved by bringing In 
the rest of the determinant. This Iteration will, no doubt, turn out 
to be co 2 n:ect If one has started out by assuming the possibilities of 
motion as the fundamental oscillation mode. However, If one has 
started out with a higher-harmonic oscillation. It will reduce to a 
suborder harmonic. How can that be prevented? 

Borkmann . - As already discussed in the research report EB 1338, 
no proof of convergence for the Iteration method mentioned exists so 
far. Just as for sjiy flutter calculation of an airplane, no more 
than three or four are singled out as decisive from the large number 
of degrees of freedom; here also a partial determinant Is singled out 
from the entire deter^nant In the expectation that the flutter mode 
of the total system will not greatly differ from that of this partial 
system. As far as this expectation comes true, the Iteration method 
offers a good chance of Improvement. Otherwise, of course, the 
application of this Iteration method may encounter difficulties. 

Lelss.- Telchmann's remarks on the sensitivity of a system to 
fluctuations In construction parameters purport that It Is not 
stiff Iclent to Investigate a system for one particular combination 
of construction parameters only but that, rather. It Is necessary 
to Include a certain variation range of the separate construction 
parameters In order to find possibly existing sensitivities. Obvi- 
ously, such sensitivities are to a large extent avoidable, but to 
that end they must first be determined. 


Translated by Mary L. Mahler 
National Advisory Committee 
for Aeronautics 
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• ziatural circular fraquancy of tha control-aurlaca toralon 
0^ - natural circular fraquency of the nlng toralon 
"achl " circular frequency of tha wing VwiyUng 



Figure 1.- Harmful and harmless frequency region (systems with three degrees 
of freedom: wing bending, wing torsion, control -surface torsion). 



0 — 


Vjjj. - critical valodty 

“ position of the auxiliary control-aurlaca canter of gravity. 

Compare sketch of system 
“schl “ circular frequency of wing healing 

- natural circular frequency of control-surface torsion 

- natural circular frequency of auxiliary control-surface 
■“ torsion 

Figure 2.- Rules for design fail to work (system in the sense of the two- 
dimensional problem with three degrees of freedom; wing bending 
control-surface torsion, auxiliary control-surface torsion) . ’ 
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schl 


critical velocity 

position of control “Surface center of gravity 
natural qlrcular frequency of tbs control-surface 
natural circular frequency of the vlng torsion 
- natural circular frequency of tbs wing bending 


torsion 


with the 
remaining 
degrees of 
freedom 
eliminated 


Figure 3 . - Variation of a construction parameter (sTStem with three degrees 
of freedom: wing bending, wing torsion, control-surface torsloi^. 
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Figure 4 .- Formula for the flutter process 
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Motlaa formula: 
Flutter equations: 

UnloownB: 
Coefficients: . 


8or(x,z;t} - ■ Fj(x,z) • el“‘ 

+ #2 ' * ••••■'■ ®n ■ Bln('''“^- " 

ai ■ + fi2 ■ B22 (v.“) + ■ 82 n(v,u) - 0, 


% • * *a ‘ Bnn(^>“> “ ° 

- amplitudes of Fj^(x,20 
Bj^(v,u) - Jpjj(v,u; o...) • Fj • dV 

- • F, • dV + /pvT(«;c...) • F, • dV + ] fpv^Cep...) • F. • dV 

vsi — ' -?v— ' -L-= -V-— i ' 


ik 


+U“(1 


ik 




*■ /Pk^"!®—) ' ^1 * ‘ F], • dO + ijp^ (v,M‘^...) ■ F^ 


dO 




v^P • ~ 


’ ^lk(v) 

Figure 5.- Flutter equations. 
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Figure 6.- Aerodynamic control -surface balance (interpretations). 
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^Uc{v) " ••••'*’ ’ T*|u^(z^dz + .... 

„ (Z) 

r' ^ V 2 

SimpUJClcatlon jfor DjjJ[z) > 

h (z) • T'jw^(z)jdz S T’(u^*) •J'Djjj{z)dz 

u * = H ^ with ^ = 
r V 2 ^ 

Assumption used as basis: 

T’^^(z)]dz Z 8q + a^Uy(z) 




z = coordinate In direction of span width 
t(z) = wing chord at point z 


Figure 7.- Simplified determinatilQn of the air-force characterlstlCB. 
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CalciilatlOQ, Uznee 




Figure 9,“ CaXculatton. tinseB rstjuljred fbr solnHoji 
of tbe flutter equations with use of the complete 
air-force formulatioiis (calculation ttmes of one 
single calculator in months) . 


Calculation times 


r 





Figure 10.' CalculaUon times required for solution 
of the flutter equations with use of simptUOed 
formulations for the air forces (calculation times 
of OnB SlligjjB CalClilBfOr glOnths) . 
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Figure 11 .- Calculation times required for solution of the flutter equations 
without limit rule for the air forces (calculation times of one single 
calculator in months) . 
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